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We apply the sea-boson method to compute the momentum distribution of a spinless 
continuum Fermi gas in two space dimensions with short-range repulsive interactions. We 
find that the ground state of the system is a Landau Fermi liquid( < Zf < 1 ). We 
also apply this method to study the one-dimensional system when the interactions are long- 
ranged gauge interactions. We map the Wigner crystal phase of this system. 

I. INTRODUCTION 

In this article we apply the sea-boson method that is now powerful enough to yield most of the well- 
known results in one-dimension, to solve for the momentum distribution of electorns in the case when the 
electrons are in a continuum in two space dimensions with short range repulsive interactions. In other words, 
the two dimensional analog of the spinless Luttinger model. As a more nontrivial application we calculate 
the momentum distribution of a spinless Fermi system in one dimension with long-range confining (gauge) 
interactions and show that the system is a Wigner crystal. The momentum distribution of this system 
exhibits some unusual features that are probably new. 

II. THE HAMILTONIAN 

Here we compute the momentum distribution of the spinless Fermi gas with short-range repulsion in two 
space dimensions in an effort to ascertain whether or not this system describes a Landau Fermi liquid. This 
exercise is simple and is an alternative to studying the Hubbard model in 2d where the algebra is quite 
involved. We expect the answers to qualitative questions such as the validity of Fermi liquid theory to 
be the same in both models since both involve very short-range interactions. We have argued before that 
whenever Fermi liquid theory breaks down, it does so maximally. That is, it breaks down for all values of the 
coupling. Thus it is sufficient to study the weakly coupled 2D system with short-range interactions where 
the analysis is straightforward. The intuition gained from this can then be transplanted to the 2D Hubbard 
model at large U, which is likely to be hard to solve even using sea-bosons. On the other hand, there may 
some qualitatively new physics in the case of electrons with spin. This will become clearer when we actually 
decide to study the Hubbard model directly in future publications. In the present case, the hamiltonian in 
the sea-boson language is given by, 

H = E C£) 4(q)^(q) + E^(-q) + 4(q)P k ' (q) + 4< (-<0] W 

kq ^ ' q#0 kk' 



1 



This hamiltonian describes a self-interacting Fermi gas provided we assume that w q = ^ 8 (A — |q|) where 
A <C kp here < v <C 1 is a dimensionless parameter. We may solve for the boson occupation number as 
follows. 



4(qM k (q)>4E7f^^ 2 (- q ) 



where A k (-q) = n F (k- q/2)(l - n F (k + q/2)). 



( \ , «qV n F (k + q/2) -n F (k-q/2) 
e RPA (q, u) = 1 + — ^ — ^ 



(2) 



(3) 



(4) 



As argued in earlier works, we have to interpret sum over modes with care so as to not lose the particle-hole 
mode, the collective mode being obvious. The is particularly important in two dimensions where we expect 
both to be present. Thus the sum over modes is defined as follows. 

duj W(q,u) /(q,w) 



J2 /(<l' w 

i 

Here the weight function is given by, 

W(q, uj) = —Im 



Jo du W(q,u) 
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(5) 



(6) 



,e_RPA(q, u - «0+) 

In our earlier work, we had suggested that in the above formula we have to use a dielectric function that 
is sensitive to significant qualitative changes in one-particle properties. The simple RPA dielectric function 
does not possess qualities we expect from a Wigner crystal. Thus we shall have to derive a new dielectric 
function using the localised basis rather than the plane-wave basis. In the case of the Fermi gas in two 
dimensions we find that the system is a Landau Fermi liquid and there is no need to use a better dielectric 
function, the simple RPA suffices. The momentum distirbution is always given by, 



< n k > = 



1 



1 + e 



-2S° (k) 



n F (k) + 



1 



-2S°(k) 



(i-Mk)) 



^(k) = ^(4_ q/2 (q)A k _ q/2 (q) 



(7) 



(8) 



S°(k)=£(4 +q/2 (q)A k+q/2 (q; 



(9) 



The computation of the boson occupation number {^Aj c (q) J 4 k (q)^ is the key to evaluating one-particle 
properties. 
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III. THE COMPUTATIONS 



In this section, we compute the various quantities defined in the previous sections. In the case of a of a 
Fermi gas in two dimensions with short-range repulsive interactions we may use the simple RPA dielectric 
function. The integrals are somewhat complicated since in two dimensions, the angular parts are very 
troublesome unlike in three dimensions. Therefore we take the easy way out and copy the results first 
derived by Stern 1 . 



, \ 1 , rnk F v q q 
eW«,«) = l + -2^-{^-C- 



( Jj_ mw\ _ l 

\2k F k F q 



-c+ 



\2k F k F q J 



} (io) 
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(13) 



g 2 (-q,w) 
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VZfi) Jo 



du> W(q,u) 
£rpa(<1iu) 



Ak(-q) • 
(- + %?) 2 



.9 2 (-q,w) 



W(q,u>) = - 2 



/■OC 

Z(q) = / duj W(q,u) 
Jo 



In general we have, 



S^(k) = 



d<7 <7 



W(g,w) f A (k,q,uj) g 2 (-q,u>) 



(2tt)2 7 - - Z{q) J a 

a r°° i z" 00 

r 

f A (k,q,iv) = 

Jo 



6(k % -k 2 -q 2 + 2kqcos(0)) 
{u) + ^ cos{e) _ eq) 2 



d6 



\2k F k F q) S 
(14) 



(15) 
(16) 

(17) 
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f M k 2 +q 2 -k 2 F + 2kqco S (0)) 

f B (k,q,u)= d9 — — — 

(w + ^cos(6») + e g ) 2 



Define, 



i(x; A,B) = j d.r 



From Mathematica™ we find, 



(i + BCos[3;]) 2 

(A-B)Tan[l 



ArcTanh 



u(x;A,B) = -2A- 



VB*-A* 



(A 2 - B 2 )V B 2 - A 2 



+ 



B Sin[x] 



(B 2 - A 2 )(A + B Cos[x}) 



Integrating by parts we find, 



f A (k 



,-2tt 

,Q,v)= / 
Jo 



h 2 - k 2 - n 2 
v 2kq 



+ cos{6)) u(6; uj - e q , — ) 
m 



,k 2 + q 2 -k 2 



f B (k,q,uj) = d9 sin9 S(- 

Jo ■zkq 



+ cos{9)) u{8;uj + e qi ^-) 



This may be rewritten as, 



k 2 - k 2 - q 2 

f A (k,q,uj) = j d9 sin6 5( F ^ 



cos(e)) u(9; uj — e q , — ) 



/ d9 sine 6{ kp * 9 -cos{6)) u{6 + -k; uj - e q ,'^) 



P7\ 
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2kq m 



e = arccos 



r fc 2 + q 2 -k 2 F 



2kq 



-k 2 -q 2 + k 2 F 



2kq 



f A (k, q, uj) = u(e ; uj - e q , — ) 

m 
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2kq 
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hp h q , 



^ fcj- fc 2 g 2 . 

Wq ' 



f B (k,q,w) = u(6 ;lj + e q , ^) 



g^ fc 2 + g 2 kp ^ ^ Q^k 2 + q 2 kp 



2kq 
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fc 2 + q 2 kp ^ ^ ^ k 2 + q 2 kp ^ 



(30) 



2fcg ' v 2/cg 

The collective mode occurs when /m[e] = 0, that is, for small enough q. This means that we have to treat 
this separately. 



1 f°° 1 f°° 



+ 



dq q f A {k,q,u c ) g 2 (-q,Lu c ) 



1 f°° 1 f 00 

5 ^ (k) = W I dqq W) Jo ^ W{q,Uj) f^ k ^9\-^) 



(31) 



+ 



(2tt) 2 



dq q f B (k,q,us c )g 2 {-(\,ui c ) 



(32) 



Here it is implicit that in W we assume that Im[e] ^ 0. The dispersion of the collective mode may be found 



using Mathematica™ . It is given below 



^ciq) = q (2tt + mi), 



(7T 2 g 2 + mirq v q + kpm 2 v 2 ) ? 



(33) 



2m 2 y/TTV q ^TT + mV q 

This dispersion is real and positive for all q and for all v q > 0. Thus in the small q limit, where using just 
the RPA dielectric function is justified and is also the limit where the close to the Fermi surface features of 
the momentum distribution is given exactly, we are justified in retaining only the coherent part. Thus we 
may write, 



1 f°° 

S A {^) ~ dq q f A {k,q,uj c ) g 2 (-q, 



1 f°° 

<$M k ) ~ J^j2 J dq q f B (k,q,uj c )g 2 (-q 7 u; c ) 

To determine whether or not Fermi liquid theory breaks down, we have to compute, 

1 f°° 

S A( k F) ~ j d 1 1 fA(kF,q,u c ) g (-q,w c ) 

1 f°° 

S%(k F )mj^^ J dq q f B (k F ,q,Lj c )g 2 (-(i,uj c ) 



(34) 



(35) 



(36) 



(37) 



If S A (kp), S B (kp) < oo then the ground state is a Landau Fermi liquid. If S A (kp) = S B (kp) = oo then the 
system is a non-Fermi liquid. For small q if we set u c — v e ff q we have, 
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f A (k F ,q,ui c ) ~ f B (k F ,q,uj c ) ~ 1/q 2 



(38) 



Also, 

g 2 (-q,u; c )~q (39) 

Thus the integrals in Eq.( 36) and Eq.( 37) are infrared finite. This means that S A (k F ), Sg(k F ) < oo and 
the system is a Landau Fermi liquid. The details of the momentum distribution can be worked out but are 
not terribly important. 



IV. ONE DIMENSIONAL SYSTEM WITH LONG-RANGE INTERACTIONS 

In this case, we expect the system to be a Wigner crystal. Thus we have to be careful about the choice of 
the dielectric function. First, we postulate that v q = 2e 2 /(q a) 2 which corresponds to the gauge potential. 
Here a has dimensions of length and e 2 > is dimcnsionlcss. From the form of this potential, one hopes 
that we need not concern ourselves with the issues that were relevant in the case of the Calogero-Sutherland 
model namely the repulsion attraction duality (more prosaically called back-scattering) . Thus we may write 
as before, 



e r2 {q,uj) + e t2 (q,u) 



/'OC 

Z(q) = / dto W(q,uj) (41) 
Jo 

9-\-*») = ~<w) (42) 

As mentioned before, we have to be extra careful in making sure that we choose the right dielectric function. 
The RPA-dielectric function is not likely to suffice since its static structure factor (SSF) does not exhibit 
the features we expect from a Wigner crystal. In particular, we expect S(2kp) = oo as we shall see soon. 
To convince ourselves of this we ascertain the properties of the RPA dielectric function with long-range 
interactions. 

2" 

(43) 



TTl 



(k F + q/2) 2 - (jfj 



( kF -q/2) 2 -(^y 

The zero of the above dielectric function gives us the dispersion of the collective modes. 
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m 



(fc F + g /2)2-(fc F - g /2)2 ea; p(-|^) 



■\ l-e*p(-^) 

For |g| <C fci? and w 9 = 2e 2 /(aq) 2 we find, 



(44) 



/ — n rrr „2l / e 2 k F m nr„2 

w c (g)w— V — \/- + " 2 2 +0(0 ) 45 

m V or V 7r ze^ra^ 

This plasmon-like gap w = m\ [^W ^\[^ m ^ e collective mode is present due to the characteristic 1/q 2 
nature of the potential. But this is also present in the three dimensional electron gas and is not a sign 
of an insulator since the latter is not at high densities. A gap in the one-particle Green function at the 
Fermi momentum could be taken as a sign of insulating behaviour 3 . However, in our approach we are 
unable to compute the full Green function as yet. Thus we must resort to a more indirect approach. For a 
Wigner crystal, the SSF must exhibit certain singularities. Thus we have to use the generalised-RPA that is 
sensitive to qualitative changes in single-particle properties. The new dielectric function will involve the full 
momentum distribution which has to be determined self-consistently using the above sea-boson equations. 
In our earlier work we suggested that the new dielectric function should also involve fluctutations in the 
momentum distribution, however it now appears that that is fortunately not needed. The number-number 
correlation function is vanishingly small in the thermodynamic limit as shown in another preprint and this 
means we may simply write, 

e(q,c) = l + ^ (46) 

and the momentum distribution is determined self-consistently using the sea-boson equations (Eq.( 7)). This 
is too difficult to solve analytically and hence we have to resort to a numerical solution. In order to simplify 
proceedings even further, we use only the collective mode. The particle-hole mode which is due to a nonzero 
Im[e] is needed if one is interested in features of the momentum distribution away from the Fermi surface 
more accurately. However we shall hope that this is given not too badly even at these regions far from the 
Fermi points. 

4(q)A k (q)) = I J^%T 2 3 2 (-q,-c) (48) 



g 2 (-q,^ c ) = 



irq 



uj 2 - {v F q + e q ) 2 us 2 - (v F q - e q ) 2 



(49) 



S°A(k) = |E , " F f*:, g) » 9 2 (-^ c ) (50) 

L q Mg) + -i - e 9 ) 2 
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cO/m 1 V" (l-n F (k + q)) 2 



(51) 



To proceed further, we have to ascertain the nature of the collective modes w c . If we use the RPA-dielectric 
function, we find a constant dispersion (plasmon) for small \q\. However we have found that this choice is 
inconsistent since if we use the momentum distribution obtained from this to solve for the dielectric function 
and recompute the collective mode we obtain a completely different answer namely : Lo c {q) = w s |q , |- Therefore 
it is critical that we get the dispersion right. It appears then that we have to use the form given in the 
appendix which is not easy to simplify. A systematic approach for obtaining the dispersion of the collective 
modes has been suggested by Sen and Baskaran 5 . Since the plane- wave basis is not appropriate for deriving 
a formula for the dielectric function of a Wigner crystal, we shall follow this approach. First, we would 
like ascertain the lattice structure in the small a limit. In this limit, the potential energy dominates over 
the kinetic energy. If we assume that the electrons are all on a circle of perimeter L then to minimise the 
potential energy, we have to maximize the separation. This leads to an equally spaced set of lattice points 
with lattice constant l c such that N l c — L. Thus we have l c = 1/po = it/hp. Thus we assume that the 
electrons all lie on a circle with equal spacing between them. Therefore we expect the structure factor to 
diverge for a momentum qo = 2ir/l c — 2kp. From the Bijl-Feynman formula S(q) = e q /uj c (q) we may suspect 
that a choice of ui c that vanishes at q = 2kp is needed. The form of the dispersion is given in the appendix. 
For tt/N <C \q l c \ <C 2tt it seems that uj q w uiq. For \q l c \ <C tt/N we have to be more careful. And of 
course we must have u> q = for q l c = ±2tt. But since in the thermodynamic limit tt/N w we may choose 
(hopefully) uj q luq. In Fig. 1 and 2 we see the momentum distribution obtained from these formulas has 
been plotted. In fact, we may write down a closed formula for the momentum distribution. 



n k = i(l + Exp 



m ljq 



h 2 — h 2 



)n F {k) + ^(1-Exp 



' U2 _ p 

A, IV jp 



)(l-n F (k)) (52) 



The striking feature of this momentum distribution is that it is perfectly flat at = kp. In other words, 
not only is the slope zero but all the derivatives of the momentum distribution vanish at |fc| — kp. This is a 
striking prediction. This may be contrasted with the smooth Gaussian function of Gori-Giorgi and Ziesche 2 ( 
Eq.(Bl) in their Appendix B ). But they consider three dimensional systems which may be different from the 
one studied here. One particle spectral functions are accessible to tunneling experiments or angle-resolved 
photoemission spectroscopy (ARPES). A more difficult problem may be to experimentally realise a Id electron 
system with long-range gauge interactions. 
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FIG. 1. Momentum Distribution of a Wigner Crystal (with zoom) 




FIG. 2. Momentum Distribution of a Wigner Crystal (no zoom) 



The formula below for the static structure factor is derived in the appendix. 

S(q) 



1 S j n 2(g^vic) 



e m "o 



This may be further simplified in the thermodynamic limit as follows. Consider, 

sin{N x) 



Then we may write, 



5(x) 



2tt 



7T X 



27T 



S(q) = 6(g)— + S(q - 2k F )— 



(53) 



(54) 



(55) 



As we can see here the strucutre factor diverges at \q\ — 2kp which means the system is a Wigner crystal ( 
a true Wigner crystal, since the divergence is from a delta-function ) 
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V. APPENDIX 



Here we use the approach suggested by Sen 5 to derive a formula for the collective modes. The formula 
Eq.( 46) although probably right is not very illuminating, for it is hard to see how the structure factor derived 
from this formula possesses the features we expect namely a divergence at \q\ = 2Uf- Thus we would like to 
derive a formula for the dielectric function where this feature is manifest. To do this we adopt the localised 
basis rather than the plane- wave basis. In real space, the hamiltonian we are studying is written as follows. 



N-l 2 

H _ >ip Pi ly 



Pi ' ' 

i=0 i>j 



(56) 



We assume that particles are on a circle and |ar| is the chord length. We would like to compute the dielectric 
function using this model. The density operator in momentum space is, 

N-l 
m=0 

Here x m = m l c + x m is measured along the circumference of the circle (see Fig. 3 below) . 



small deviations 




lattice point 



FIG. 3. Schematic Diagram of Electrons on a Circular Lattice 



The average density is given by, 



JV-l 



(p q > = eiq m lce ~ iq 



(58) 



m— 



It can be shown that (see below) < >= l/me^o, independent of the index n. Thus we have, 



(p q > = - 



1 _ e *gN l c 



oiq lc 



(59) 



In those instances where < p q >^ the static structure factor is given by, 



S(q) = ' < Pq > 



N 



N 



I _ e iqN l c 



1 - e l i 



in l c 



2_ sin 2 (^) 
N S m 2 (SA) 



e 



(60) 
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The rest of the details are as follows. We write x m (t) w R m (t). In terms of the small angles 9i the 
hamiltonian in Eq.( 56) may be written as follows. 



H = 



1 ^ 8 2 Re 2 



2m R 2 ^ mi 2a 2 ^ 



( cos ( 27r -+0 m )- cos (2^— +0 m O) 2 + (^(27r-- 1 - ■'- - 



l )_ ain( 2 7r _+0 m ,)) 2 



(61) 



We may expand the above hamiltonian in powers of the angle and retain only the leading terms to arrive at 
the following hamiltonian in the harmonic approximation. 



JV-l _ 2 



where, 



A(n, n ) = 



7re 

4ia 2 



■ , (n-n) 
— ]y — ) 



(62) 



(63) 



z — 71 1 f71 — ft ) 

B(n,n ) = «sn(**»(f — — — )) cos(n — — ) 



(64) 



Despite appearences to the contrary, the extremum of the potential is at x n = 0. Since A > 0, this extremum 
is also a minimum. One has to now compute the various correlation function of the system. The primary 
one of interest is, 



Gu(nt; nt)=<T x n {t) x n > (t ) > 



(65) 



The other is, 



G 2 i(nt;nt) =<Tp n (t) x n >{t ) > 



(66) 



Thus we have, 



d i 
i—G\\(nt; n t ) = — G 2 \(nt\ n t ) 
at m 



(67) 



i^-G 2 i(nt- n't) = 6 nn ,S(t - t) - 4i ^ A(n,j)(G n (nt; n't) - G n (jt; n't')) 

This may be solved by a Fourier transform. 

G^nt-nt) = * £yp('-'')l£y« '-("-^Gyfozp) 
v <? 



(68) 



(69) 



Thus we have, 



i z p G u (q,Zp) = —G 21 (q,z p ) 



m 



(70) 
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Thus, 



i z p G 21 (q,z p ) = l + 4i (A(q)-A(0)) Gu(q,z p ) (71) 

%) = E% i,iei (72) 
j 

Gn («, z P ) = (i m z 2 p - Ai (A(q) - i(0))) (73) 

Gu(nt;n't) = ^E e<<? Wn ""' ) jE e * p( *" t ' ) ( TO *p " 4 " A (°») ( 74 ) 

p 



JV-l 



A(a) = £° V fe^+« '«W - e J (-*+9 'cW^ = ££ f 1 + gf ^ 1 ± f!! — ) (75) 

i(0) = ? f— - 2 - ^ (76) 

w 2iVl-e 4 ^ 1-e-^y 27ra 2 v ; 

Here Co = 7re 2 /4ia 2 . The dispersion of the collective mode is then given by, 

^ = \[^{M0)-Mq)Y (77) 

If q = or q = 2kp then oj q = 0. For a more thorough analysis one has to compute the full dielectric function 
from the DDDCF and use it to compute the full momentum distribution that is accurate even away from the 
Fermi surface. However we shall be content at features close to the Fermi surface. For the static structure 
factor we have to compute the equal time version of the correlation function. We find that iv q is in general, 
complex. This means that the eigenmodes also have a finite lifetime. Thus we have, 

■t , -1 />oo 

Gll(nt;n 'i) = _£ e * '=(«-»)__ jj Zp (zl+^q))- 1 



1 p iq l c (n-n ) 



2m N ^ Lu c (q) 

From Eq.( 77) it is clear that for ir/N <C \q l c \ < 2ir we have w q w ujo since A(q) « for q in this region. 
Since in the thermodynamic limit, this is all of q, we shall boldy write, 
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